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ABSTRACT 

Wc present measurements of the cosmic microwave background (CMB) power spectrum made by 
the Atacama Cosmology Telescope at 148 GHz and 218 GHz, as well as the cross-frequency spectrum 
between the two channels. Our results clearly show the second through the seventh acoustic peaks in 
the CMB power spectrum. The measurements of these higher-order peaks provide an additional test 
of the ACDM cosmological model. At £ > 3000, we detect power in excess of the primary anisotropy 
spectrum of the CMB. At lower multipoles 500 < £ < 3000, we find evidence for gravitational lensing 
of the CMB in the power spectrum at the 2.8cr level. We also detect a low level of Galactic dust in 
our maps, which demonstrates that we can recover known faint, diffuse signals. 
Subject headings: cosmology: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 

Accurate measurements of the arcminute scale tem- 
perature anisotropies in the mm-wave sk y are uncov- 
ering a complex, yet r evealing picture ([Luckcr et al.l 
12010; iFowler et all I2010L hereafter F10). On interme- 
diate scales (500 < I < 3000), the primordial acous- 
tic features imprinted on the cosmic microwave back- 
ground (CMB) at last scattering (z ~ 1100) dominate, 
with subtle distortions expected from gravitational lens- 
ing by intervening large-scale structure. This interme- 
diate range of multipoles is often called the damping 
tail of the CMB, as the acoustic oscillations are expo- 
nentially damped due t o photon diffusion (jSilkl 119681 : 
iBond fc Efstathioul[l98i ). On the smallest scales [l > 
3000), emission from radio and dust-enshrouded star- 
forming gal axies, together with the Su nyaev- Zcl'dovich 
(SZ) effect (jZeldovich fc Sunvaevl [i960 ) dominates over 
primary CMB fluctuations. 

The damping tail measurements are an additional test 
of the predictions of the ACDM cosmological model, 
a model t hat is an excellent fit to current CMB 
data (e.g.. lLarson et all 12010 : iRcichardt et alJ l2009bb 
iBrown et al.l 12000). large-scale structure measurements 
(e.g.. iReid et alJ |2010V iPcrcival ct al.l 120101). supernova 



obser vations (e.g.. iKesslcr ct al.l I200S : lAmanullah et al.l 
l2010f l , an d a host of other astro nomical observations 
(see, e.g., iSpergel et al.1 l200~7l and iKomatsu ct all 120101 

for reviews). The amplitude of the fluctuations in the 
damping tail is a sensitive probe of matter fluctuations 
at k ~ 0.1 — 0.25 Mpc -1 — thus, precision measure- 
ments constrain the spectral index of the primordial cur- 
vature perturbations, n Sl and its variation with scale. 
Because the positions of the high order acoustic peaks 
are sensitive to the evolution of the sound speed of the 
universe and its composition, these measurements con- 
strain the primordial helium fraction and the number 
of relativistic species including neutrinos (see e.g. IWhitd 
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Ichikawa ct al.ll2008aU bl; IKomatsu et al.ll2010ft . While the 
damping tail gives us leverage on early universe cosmol- 
ogy, the composite tail (£ > 3000) is sensitive to a variety 
of astrophysical phenomena, including the properties of 
the intracluster medium, the redshift distribution an d 
clustering of dusty sub- mm gala xies (lHall ct al.l I2010D 
and t he physics of rcio ni zation (iHuffcnbcrgcr Sz Seliakl 
[20051 ISehgalcna1l27)07l ). ISehgal ct al.l (I2010bl ) describe 
theoretical expectations for these small-scale measure- 
ments. 

In this paper, we present a measurement of the 
CMB power spectrum over the range of multipoles 
500 < i < 10000 from the Atacama Cosmology Tele- 
scope (ACT) using the 148 GHz and 218 GHz chan- 
nel data from the 2008 observing season . ACT is a 
mm-wave, arcminute- res olution telescope (jFowlcr et al.1 
l2007t iSwctz ct al.l 120101 ) custom built to make pre- 
cise observations of the microwave sky over the damp- 
ing and composite tail regimes. In recent years, sev- 
eral groups have reported rapidly improving measure- 
ments of fluctuations over various portions of this 
multipolc ran ge : Bol ocam dSavcrs et al.ll2009D . QUaD 
(IBrown et all 12000; [Friedman et al.l I2009D. APEX- 
SZ (IRcichardt ct all l2009al). AC BAR (IRcichardt et al 



20061 ) CBI (ISicvcrs et al.l [2009h . S PT dLuckcr ct all 
27)10 lHall ct al.H20lO) and ACT (jFlOl ). The current pa- 
per enhances the results of iFlO in a few ways. First, 
we augment our 148 GHz data with results from our 
218 GHz channel. Second, the power spectrum estima- 
tion methods have been revised to provide enough angu- 
lar frequency resolution to detect the acoustic features 
on the damping tail. Third, the area used for the power 
spectrum analysis has been increased from ~220 deg 2 to 
~ 300 deg 2 . 

This paper is one in a set of papers describing ACT 
and its 200 8 Southern survey. The ACT instrument is 
described in lSwetz ct alJ (I2010D. a 148 GHz p oint source 
catalog is presented in iMarriagc ct all (|2010bD (hereafter 
M10), and a 148 GHz S Z cluster catalog is pres ented in 
iMarriage et all (|2010af ). iMcnanteau ct al.1 ([2010 ) discuss 
the multi-wavelengt h followup of ACT clusters, while 
ISehgal ct all &)]M i present as constraints from SZ clus- 
ter de tections. On the power spectrum side. lHaiian et ahl 
(|20l0 ) report on the calibration of ACT maps using 
cross-correlations wi th WMAP seven-year maps, and 
iDunklcv ct al.l ([2010) present the constraints on cosmo- 
logical parameters derived from the power spectrum pre- 
sented here. 

The paper is organized as follows. In Section [21 wc 
briefly review the ACT instrument and observations 
made so far, touching upon mapmaking and beam de- 
termination techniques, and the calibration of our data. 
In Section [3) wc describe the method used in power spec- 
trum estimation. Simulations used to test our pipeline 
are discussed in Section 2J We present our results in Sec- 
tion[5]and discuss tests performed to validate our results. 
The main sources of foreground contamination and the 
methods used to treat them are discussed in Section [6j 
We discuss the lensing contribution to the power spec- 
trum in Section before concluding in Section [8] 

2. OBSERVATIONS AND MAPMAKING 

ACT is a 6-m eter off-axis Gregorian telescope 
(| Fowler et al.ll200~7l ) situated at an elevation of 5190 me- 
ters on Cerro Toco in the Atacama desert in north- 
ern Chile. ACT has three frequency bands centered 
at 148GHz (2.0mm), 218GHz (1.4mm) and 277GHz 
(1.1 mm) with angular resolutions of roughly 1.'4, l.'O and 
0.9, respectively. The high altitude site in the arid desert 
is excellent for microwave observations due to low precip- 
itablc water vapor and stability of the atmosphere. The 
tropical location of ACT permits observations on both 
the northern and southern celestial hemisphe res. Further 
detail s on th e instrument are presented in ISwctz et al.l 
(|20l0 ). IFlO and references therein^. As of this writ- 
ing, ACT has completed three observing seasons (2007, 
2008, 2009) surveying two stripes on the sky: a 9°-wide 
stripe centered on declination -53?5 and a 5°-wide stripe 
centered on the celestial equator. The power spectrum 
presented here is derived from the southern stripe data 
from the two lower frequency channels in 2008. An im- 
portant aspect of ACT's scanning strategy is the cross- 
linking of observations. Every point in the survey area is 
scanned along two different directions during the night. 
We use constant elevation scans centered on two different 
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azimuths — once when the survey area is rising on the 
sky, and again when it is setting toward the end of the 
observing night. In principle, this cross-linking allows 
an unbiased reconstruction of all modes in the map, and 
permits recovery of the power spectrum for multipoles 
as low as a few hundred with errors completely domi- 
nated by cosmic variance. A more detailed a ccoun t of 
the o bserving strategy of ACT can be found in iFld and 

ImIoL 

2.1. Mapmaking 

For the 148 GHz data, we use the maps from IF Id . in 
which details of the data reduction and mapping can 
be found. We model the data as d = Pm + n, where 
P projects the map m into the time stream, and n is 
the noise, with covariance matrix N. There are no con- 
straints on the contents of P and m, which can contain 
multiple components. For ACT, in addition to solving 
for the map of the sky, we solve simultaneously for noise 
correlated amongst the detectors. Under this model, the 
data can be described by 

d = Pm + Ac + n (1) 

where A are (assumed constant) patterns of correlation 
across the array (such as a common mode) and c are the 
timestreams associated with each pattern in A. This is 
mathematically equivalent to having separate blocks in 
a generalized projection matrix and a generalized map 
solution. The maximum likelihood solution, m, is then 
given by solving the standard mapping equation: 

P T N _1 Pm = P T N _1 d. (2) 

We solve for in iteratively using a preconditioned 
conjugate gradient (PCG) method (|Press "eTall [2003 ; 
iHinshaw et al]|2007f >. 

The higher atmospheric noise in the 218 GHz data re- 
quired some changes to the mapping pipeline. For the 
148 GHz data, we found that taking the array patterns 
A to be the eigenvectors corresponding to the 10 largest 
eigenvalues of the data covariance matrix for each 15- 
minute chunk of time-ordered data (TOD) worked well 
for correlated noise rejection. For the 218 GHz data, 
substantial atmospheric power remains with this pre- 
scription, and so we altered it as follows: first, for each 
TOD we take the band-limited data between 0.25 and 
4 Hz, find the eigenvalues and eigenvectors of the cor- 
responding covariance matrix, and keep all modes with 
eigenvalues larger than 3.5 2 times the median eigenvalue 
(equivalently modes with timestream amplitudes larger 
than 3.5 times the median). We typically find 30-50 
modes. Then, we create the covariance matrix from the 
data high-pass filtered above 4 Hz, project out the modes 
already found in the 0.25-4 Hz band, and keep all re- 
maining modes with eigenvalue larger than 2.5 2 times 
the median. We typically find 1-2 additional modes in 
this step. Of the several different mode removal schemes 
we tried, we found this fairly aggressive one gave the best 
signal-to-noise on intermediate and small angular scales, 
where the 218 GHz data are most valuable, at the price of 
worse signal-to-noise and slower convergence of the map- 
per on large scales. Since the method estimates both the 
correlated modes and the map of the sky simultaneously, 
mode removal does not bias the maps, although it makes 
some sky map modes noisier. 



The second change in the mapping is the use of a 
preconditioner with the conjugate gradient technique. 
In general, a preconditioner is an approximation to 

(P T N _1 P) that is used to speed the convergence of the 
conjugate gradient solution. In general, we would use a 
two-piece preconditioner, with separate parts for the sky 
solution and the correlated noise component. However, 
if the columns of A are diagonal (which they are by con- 
struction), the natural preconditioner for the correlated 
part of the noise is simply the identity matrix and so we 
neglect it. The sky preconditioner is a diagonal matrix 
with elements equal to the inverse of the "n bs" map — 
a map of the number of observations in each pixel. We 
find the use of this preconditioner greatly aids in the con- 
vergence of the conjugate gradient solution for 218 GHz 
data, unlike for the 148 GHz data where it makes little 
difference. 

2.2. Beam Measurements 

The window function due to the beam shape is 
obtained in a m ethod similar to that followed in 
iHincks et all 120091 hereafter H09. We use 22 (26) maps 
of Saturn from throughout the 2008 season to estimate 
the beam shapes for 148 GHz (218 GHz). The maps are 
mad e with an independent pipeline from that used in 
IH09I . but produce consistent results. Each map is binned 
into a symmetrized radial beam profile, with 0.'17 bins, 
out to a radius, < ax , of 8'(6') for 148 GHz (218 GHz).' 
The dominant noise source in the maps comes from the 
slow variation of atmospheric brightness, and since this 
projects into relatively large angular scales it causes sig- 
nificant noise correlation between radial bins in each pro- 
file. The individual profiles are thus used to estimate the 
mean profile, and the full covariance matrix between bins 
is estimate d from the scatter between them. This is in 
contrast to IBM where the mean profile was computed 
from a stacked map, and only diagonal errors were con- 
sidered. The full covariance formulation serves to prop- 
agate the large scale map noise into law-£ uncertainty in 
the window function. The zero level of the radial pro- 
file is poorly constrained in the maps due to atmospheric 
contamination at large radii. Because the sharp circular 
edge of the cold Lyot stop produces an Airy pattern at 
large scales, the radial profile is expected to fall as 1/0 3 
so a fit to this form beyond estimates the 

zero level and its covariance with the binned profile. The 
resulting function is also used to extrapolate the beam 
when computing the beam transform. 

Following the same prescription as IH09I . a small set of 
basis functions is fitted to the inner part of the beam 
profile. A natural choice of basis functions in harmonic 
space are Zernike polynomials which are compact on the 
unit disk, or their transforms in angular space which are 
Bessel functions divided by a linear function of radius. 
For 148 GHz (218 GHz) only 25 (18) basis functions are 
required to approximate the beam profile at the 1% (2%) 
level within 1' and below the 15% level within 8'( 6'). 
Using a small set of basis functions simplifies calcula- 
tion of the covar iance matrix of the window function 
(jPage ct al.ll2003D from the covariance of the coefficients 
of the basis functions and the wing fit. For the pur- 
pose of calibration against WMAP maps ()Haiian ct alj 
I2010T ) . the normalization of the window function is fixed 
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Fig. 1. — The rms temperature uncertainty for one-arcminute pixels for the 148 GHz maps (top) and 218 GHz map (bottom). Also shown 
in bold lines are the boundaries of the four rectangular patches used for spectral analysis. 



to unity at £ = 700, and the calibration error is fac- 
tored out of the covariance matrix (described in H09) 
so that the window function has zero error at £ = 700. 
Between £ = 1000 and 10000, uncertainties in the win- 
dow function at 148 GHz are between 0.7 and 0.4%, and 
at 218 GHz between 1.5 and 0.7%. The uncertainty in 
the window function is significantly higher for £ < 700, 
where atmospheric contamination makes measurements 
very difficult. In this paper, we use the beam window 
function only. The full covariance matri x of the window 
function is used in iDunklev et al.l (|2010D . 

2.3. Calibration 

ACT scans every point in the survey area both when 
it is rising and setting. This cross-linking along with 
the unbiased map-making method described in Sec- 
tion 12.11 allows the reconstruction of all modes in the 
maps, wi thout biasing the la rge-angle modes. As de- 
scribed in lHaiian et al.l (|2010D . we calibrate the 148 GHz 
ACT maps directly to WMAP 7-year 94 GHz maps 
(| Jarosik et all [20Toh of the identical regions. Thus, the 
cosmological analysis is done with the same data used 
for the calibration. By matching the ACT- WMAP cross- 
spectrum to the ACT p ower spectrum and the WMAP 
7-year power spectru m (lLarson et al.ll2010D in the range 
400 < £ < 1000, lHaiian et all (|2010D calibrate the 
148 GHz ACT spectrum to 2% fractional temperature 
uncertainty. Similar methods are applied to 218 GHz 
maps, but their larger map noise levels on large angular 
scales mean that the result does not currently improve 
on the results of an independent method of calibration, 
based on observations of Uranus^, described below. 

2 It may be noted that although we use observations of Saturn 
for evaluating beam shape, we do not use them for calibration. Re- 



ACT made approximately thirty usable observation of 
the planet Uranus during its 2008 season. We gener- 
ated a map per observation after calibrating the time- 
ordered data to detector power units, and determined 
the peak response of the planet with corrections for tem- 
perature dilution due to the finite instrumental beam 
size. The result is then compared to the Uranus tem- 
perature in CMB differential units at the effective band 
cente r to obtain a calibration factor (sec lHaiian et al.l 
mm for more details). We take the brightness temper- 
ature of Uranus to be 107 16 K and 96 ± 6 K for the 
148 GHz and 218 GHz bands respectively. These temper- 
at ures are based on a repr ocessing of the data presented 
bv lGriffin k Ortonl (|1993l) . in combination with WMAP 
7-year measurements of Mars and Uranus brightnesses 
(jWciland et al.ll2010D . We find an uncertainty of 7% in 
the Uranus-based calibration that is dominated by the 
6% uncertainty in the planet's temperature. The ab- 
solute calibration is consistent with the WMAP-based 
calibration described above. For 218 GHz, we adopt the 
Uranus-based calibration. 

3. POWER SPECTRUM METHOD 

iDas ct all ()2009D describe the basis of our power spec- 
trum method. Th is paper improves on the analyses of 
iFowlcr ct al.l (|2010T ) by optimizing the angular frequency 
resolution of the spectrum over the damping tail to re- 
solve the acoustic features, as well as by extending the 
multipole range of the power spectrum to I — 10000. 

3.1. Fields used for Spectra 

liable modeling of the brightness temperature of Saturn is compli- 
cated by the effects of its rings. The effective temperature changes 
with the ring inclination over the season. This is not the case for 
Uranus, whose brightness temperature is fairly stable. 
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Fig. 2. — Beam deconvolved noise power spectra of the 148 GHz 
(open circles) and the 218 GHz (filled circles) maps shown against 
a theoretical lensed CMB spectrum (solid line). The theoretical 
white noise levels are shown using dashed (for 148 GHz) and dotted 
(for 218 GHz) lines. 

The power spectrum is computed from independent 
analysis of four contiguous patches in the ACT southern 
stripe (see Fig. [1]) that together cover a rectangular area 
of 296 deg 2 from a = 00 h 22 m to 06 h 52 m (5?5 to 103°) in 
right ascension and from 5 = —55° to —50° in declina- 
tion. Each patch is 5° x 14?8 in size. The full dataset is 
also split into four subsets by distributing data from each 
of four successive nights into different subsets, thereby 
producing four independent maps of identical coverage 
and roughly equal depth. We refer to each subset as a 
split. Therefore, for each of the two frequency bands, we 
have four patches with four splits in each patch. Each 
split is properly cross-linked, i.e., it represents the maxi- 
mum likelihood solution of a dataset containing both ris- 
ing and setting observations. A typical white noise level 
in the 148 GHz season map (all four splits taken together) 
range from 25-50 [a.K-arcmin, while for the 218 GHz map 
it varies between 40 and 70 |xK-arcmin, with the noise 
being higher toward the edges of the maps (Fig. [1]). The 
noise power spectra of the 148 GHz and 218 GHz maps 
are shown in Fig.[2j At small scales, the noise is typically 
white, while at larger scales there is a plateau of noise 
with a red tail towards low multipolcs, which is primarily 
attributable to atmospheric contamination. 

3.2. Preprocessing of Maps 

Before patches are cut from the splits, each split map 
is high-pass filtered in Fourier-space. The high-pass fil- 
ter is designed to remove all modes below £ = 100, and 
suppress the modes between £ = 100 and 500 with a sine- 
squared function, which rises from zero to unity within 
that range. The very lowest {I < 100) multipoles are 
dominated by extremely large atmospheric noise, and 
rolling off the filter to £ = 500 prevents these modes 
from leaking into and unnecessarily biasing the higher 
multipole portion of the spectrum. 

Next, we prewhiten the maps. Prewhitening is a lo- 
cal, real-space operation on the map designed to reduce 
the d ynamic range of its Fourier components (|Das et alj 
120091 ). Prewhitening is particularly important for the 
current spectrum, as it is designed to specifically target 



the damping tail of the CMB (1000 < £ < 3000). This 
region of the spectrum has a steep slope (Ce ~ £~ 4 ). 
In the absence of prewhitening, finite boundary effects 
and application of the point source masks cause a large 
amount of power to be aliased from low to high mul- 
tipolcs contaminating the spectrum at large multipolcs. 
Although an unbiased estimate of the spectrum is ob- 
tained by deconvolving the mode- coupling window from 
the spectrum (See Section l3~5l and [Das et al.ll2009l .). this 
contamination adds to the uncertainty on the estimate, 
resulting in unnecessarily large error bars. Prewhitening 
is performed by taking the difference of two versions of 
the same map, one convolved with a disc of radius 1' 
and the other with a disc of 3' . To simulate discs with a 
pixelated kernel, we use the cloud-in-cell approach where 
we assign appropriate weights for pixels under the disc 
according to its area of overlap with the disc. This op- 
eration approximately amounts to taking the Laplacian 
of the map, which effectively multiplies the map Fourier 
transform by ~ (£R) 2 and the spectrum by ~ (£R) A in 
the multipole range 900 < £ < 3000. Here R = 1' is the 
radius of the smaller disc. This flattens out the damping 
tail but makes the less colored lower multipole portion 
of the spectrum steeply rising. This transition occurs 
around £ ~ 900, where £ 2 R 2 0.02. Therefore, we add 
2% of the original map back to the resulting map and this 
flattens out the £ < 900 portion of the spectrum. Note 
that the prewhitening kernel is designed to flatten out by 
£ r-j 3500 and therefore leaves the high multipole compos- 
ite tail (which is already white) largely unaffected. 

We perform these actions in real space, but calculate 
the resulting prewhitening transfer function in Fourier- 
space so we can undo the prewhitening in the final map 
estimate. Prewhitening reduces the scatter in the £ = 
2000 — 4000 region of the power spectrum, which provides 
significant statistical weight for extracting information 
on the SZ and correlated point source signals. 

3.3. The Data Window 

After the maps are prewhitcned, four patches are cut 
from each split. Each split patch is then multiplied with a 
window. The window is a product of three components 
- a point source mask, the ?i b s map, and a tapering 
function. The details of the point source mask are given 
in Section [6] To simplify the application, we generate 
a single n Q b s map per patch for each frequency by sum- 
ming over the four splits for that patch. For the cross- 
frequency spectra we again generate a single n b s map by 
summing the ?i b s maps from the two frequencies. Mul- 
tiplication by the n Q b s map essentially downweights the 
relatively scantily observed and poorly cross-linked pixels 
toward the top and bottom edges of the map. Finally, 
to avoid ringing from the sharp edges of a patch, each 
patch is multiplied by a taper that gently rolls to zero at 
the edge of the map. We use a simple taper generated 
by taking a tophat function which is unity in the center 
and zero over ten pixels at the edges, and convolving this 
with a Gaussian kernel of FWHM ~ 2.5'. 

We will denote the two frequency channels as A and 
B, and use i,j,k,l to denote the sub-season splits. The 
patches will be denoted by Greek indices. We will dis- 
tinguish the data windows only by their patch indices, 
and denote them by W Q (n), suppressing the channel la- 
bels. Thus, split i of the patch a of frequency A, after 
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windowing, become^: 



T A {h) = W a (n)T™(n 



iA 



(3) 



where Tp5 a (n) denotes the map after filtering and 

prewhitening, and is related to the original map T^ A (n) 
through the Fourier-space relation: 



T A (i)=F i T a A {l), 



(4) 



where Fg is a product of the prewhitening and beam 
transfer functions, the pixel window function and the 
high-pass sine-squared filter. Here, and in what follows, 
the boldface i is meant to represent a 2D wave-vector 
in Fourier-space, while we use regular £ to represent its 
the absolute value. Symbols subscripted by £ represent 
quantities that are isotropic, while ones shown with an 
argument (£) or subscripted by £, stand for explicitly 2D 
quantities in Fourier-space. 

3.4. Azimuthal Fourier-mode Weights and Binning 

After windowing the data, a 2D pseudo-spectrum is 
computed as: 



iAxjB _ 



Re 



rp*iArpj 1 



(5) 



where the patch index has been suppressed to simplify 
notation. At this stage, we crop out and retain a rect- 
angular area of the 2D spectrum defined by —10000 < 
4 < 10000 and -10000 < £ y < 10000, throwing out 
all information at £ x ,£ y > 10000. This trimming re- 
duces the number np of Fourier-space pixels by a factor 
of ~ 4, and is crucial for several subsequent operations, 
especially the direct computation of the mode-coupling 
matrix (see Section [3. 5J) , which involves a computational 

step that scales as np. The ID binned spectrum C\ re- 
sults from averaging the 2D spectrum in annular bins: 



^-liAxjB _ 
°6 



iAxjB 



(6) 



where P^t is the binning matrix, whose value is zero when 
t is outside the annulus defined by bin index b. Because 
the noise properties are not isotropic in Fourier space, 
we employ a 2D weighted average, with the weights de- 
termined by the noise power spectrum. We briefly de- 
scribe here how the weights are computed. For a given 
frequency, say, A, we compute the mean 2D auto spec- 
trum from the four split patches and subtract the mean 
2D cross-power spectrum from it. This gives an esti- 
mate of the 2D noise power spectrum for the split patch 
N AA . We assume that each split has roughly equal noise, 
and approximate the season noise power spectrum as 
Nf A /4 . We estimate the 2D variance of the spectral 

estimator Cl Ax -' B as 



a 2 (Ci AxjB ) ~ (Cf bf + N AA /A )(Cf bf 2 



iVf B /4), 
(7) 



3 Throughout the paper, we use X(fi) to denote a real space 
quantity, and X(£) or X% to denote its Fourier transform. We 
use the tilde to denote windowing (not Fourier transformation), 
so that X(h) is W(n)X(ri) after multiplication by a window. All 
calculation is done under the fiat-sky approximation. 



where Cf 1 is an azimuthally symmetric theoretical CMB 
spectrum, and b[ represents the beam transfer function. 
The first estimate of the azimuthal weights is chosen to 
be the inverse of the variance (Eq. [7]): 



iAxjB 



1 



iAxjB \ 



(8) 



The resulting 2D Fourier-space map^\ of weights is typi- 
cally noisy, and before using it to bin the power spectrum, 
we treat it in the following manner. In every annulus, 
defined by the bins 6, we identify outlier pixels whose 
values are greater than a threshold a& (typically chosen 
to be 10 times the median in that annulus) and set their 
values to «&. We then divide the pixels in each annu- 
lus by their mean, and smooth the resulting 2D Fourier- 
space map with a Gaussian with FWHM of three pixels. 
These operations are designed to produce well-behaved 
Fourier-space weight maps, with emphasis on bringing 
out the overall azimuthal variation in the weights by 
removing outliers and radial dependence. The Fourier- 
space weight map used in the 148 GHz x 148 GHz spec- 
trum is shown in Fig. [3] A comparison of this Fourier- 
space weight map with the signal and noise levels in 
Fig. [2] shows the expected behavior of the weights being 
isotropic in the signal-dominated regime at low £, with 
the anisotropy becoming more pronounced as the noise 
begins to dominate. The main sources of anisotropy in 
the weights are the rays of excess noise along and per- 
pendicular to the scan directions, leading to two slightly 
rotated (because the rising and setting scans are not ex- 
actly orthogonal) "X" shaped patterns of low weight re- 
gions. Another significant noise term is attributable to 
scan synchronous noise, likely caused by instabilities in- 
duced by acceleration at scan turn-arounds. The noise 
leads to horizontal striping in the maps and manifests it- 
self as an excess of power in a vertical strip —90 < £ < 90 
in the Fourier space. We set the pixels inside this strip 
in the 2D Fourier-space weight map to zero. In terms of 
the 2D weights, the binning matrix can be expressed as, 



Phi = 



w i 



(9) 



eeb 



It is noteworthy that both the n b s weighting in real 
space and the azimuthal weighting in Fourier space lead 
to a 5-10% percent lowering of the uncertainties in the 
final spectral estimates. 

3.5. Mode-coupling 
From Equations ([3]) to we can express the binned 
pseudo-spectrum C& in terms of the underlying spectrum 
CV: 



iAxjl 



(10) 



This quantity can be related to a binned version of the 
the true spectrum C& through an inverse binning opera- 

4 The 2D Fourier transform X(£) of a real-space map X(h) is de- 
fined on a 2D pixclatcd grid in ^-space. We use the term "Fourier- 
space map" to denote any 2D quantity defined on the same £-spacc 
grid. Such a quantity can be solely defined in £-space and is not 
necessarily derived from a Fourier transform of some known real- 
space quantity. 
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Fig. 3. — Azimuthal weights used for binning the 148 X 148 GHz 
spectrum. The horizontal stripe is an artifact of the symmetry of 
the Fourier space, and does not affect our calculations. The vertical 
stripe of zero weight is added in by hand, as described in the text. 
Due to the symmetry of Fourier transforms, only the upper half 
plane of this weight map is independent. The small grainy disk 
at the center corresponds to £ < 500 where we do not perform 
any smoothing — these modes are discarded from the final power 
spectrum. 

tor Qa, which is unity when t £ b and zero otherwise: 
Ct AxjB = £ Pu\W(£-£')\ 2 F?,Q e , b ,Cl AxjB 

l,t',V 

= Y^M bb >Ci A * jB , (11) 
v 

where Mbv is the mode-coupling matrix. We compute 
the Mbf exactly without resorting to any of the one- 
dime nsional approximat i ons often used in the flat-sky 
case Jffivon et all [2001 iDas et al.l 120091: iLucker et at] 
l2010T l . The mode-coupling matrices are well behaved and 
stable to inversion. We define the unbiased, decorrelated 
estimator of the power spectrum (indicated with the cir- 
cumflex) : 

QU*iB ^Mtfd**'* . (12) 
b' 

For a single frequency spectrum, we compute the cross- 
power spectrum in each patch as the mean of the six 
cross-power spectra (from four splits). We also com- 
pute the auto-spectrum as the mean of four auto-spectra 
from the splits. The variance for a given patch is de- 
termined by combining the noise estimate (obtained as 
one-fourth the difference of the mean auto and the mean 
cross-spectrum) in that patch with a theoretical signal 
power spectrum. We combine the four mean cross-power 
spectra from the patches with inverse variance weighting 
to obtain the final spectrum. 

For the cross-frequency power spectra, we compute 
twelve cross-power spectra in each patch by crossing each 
split for a frequency with the three other splits from the 
other frequency. We do not cross the same splits (con- 
sisting of data from the same nights) between the two 
frequencies because of the possible contamination from 
covariant atmospheric noise among the two channels. We 
proceed by averaging those twelve per patch, and then 
combining the cross-spectra from the patches with in- 
verse weighting as above. In this case, we estimate the 



variance in a patch as the square root of the product of 

the single-frequency variances i n that patch. 

For parameter estimation fscc lDunklcv et aT1l2010D . we 

wish to compare our spectral estimates Cb to theoretical 
power spectra Cf 1 . To do so we find the bandpower win- 
dow function, Bbg which converts Cf v to binned theoreti- 
cal spectra Cj h . (We suppress the superscripts on the Cs 
for clarity.) We begin with Equation ([TO]) which relates 
the binned spectrum to the underlying 2D spectrum, and 
introduce a matrix Iu that describes the mapping of a 
theoretical spectrum defined at integer multipolcs onto 
our discrete Fourier pixels. The binned theoretical spec- 
trum follows from Equation (|12p : 

b',£",e' 

= Bb£ Cf 1 . (13) 

Since each rectangular Fourier-space pixel spans a range 
of integral multipoles, it gets contributions from several 
multipolcs in the theoretical spectrum. We choose la as 
a Gaussian with its peak at the pixel center (£ x ,£ v ), and 
width corresponding to the pixel size, so as to give more 
weight toward the center of the pixel. The bandpower 
window functions are insensitive to the exact form of this 
response function. As with the mode-coupling matrix, we 
compute Bu exactly. 

3.6. Summary of power spectrum method 

The angular power spectrum is conceptually simple, 
but measuring it over a wide range of angular scales 
with sufficient angular resolution to see acoustic peaks 
from maps over a small portion of sky with significant 
variation in sensitivity over the map, is technically chal- 
lenging. As a quick reference, we summarize here our 
multiple-step procedure to accomplish this task: 

• Divide of the data into independent "splits" each 
comprising a quarter of the total observation 
nights, and construct independent maps for each 
split; 

• High-pass filter the maps, to eliminate the large 
angular scale modes which are poorly constrained 
and can bleed power into the smaller-scale modes; 

• Prewhiten the maps through real-space convolu- 
tions so that the power spectrum of the maps is flat- 
tened, which reduces aliasing of power from large 
scales to small scales given the steep power spec- 
trum of the maps; 

• Divide each split map into four independent 
patches and construct a window for each patch 
which accounts for varying sensitivity over the 
patch, perform masking of point sources, and edge 
tapering to avoid spurious power from a sharp cut- 
off; 

• Calculate 2D binned pseudo-spectra from 
the Fourier transformation of the windowed, 
prewhitened, filtered maps, for each pair of 
frequency, split, and patch values; 
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Fig. 4. — Anisotropic noise in data and simulated patches. Top: Difference map of two splits in an 148 GHz data patch (left), and its 
2D power spectrum (right). Bottom: Same as above, except for a simulated patch. This random realization was seeded by the noise power 
spectrum of data patch shown in the top panel, as described in the text. 



• Construct a onc-dimensional binned pseudo- 
spectrum via azmiuthal averaging of the two- 
dimensional spectrum over bins in multipole £; 

• Estimate the true binned power spectrum from 
the binned one-dimensional pseudo-spectrum by 
inverting a mode-coupling matrix which accounts 
for the effects of beam profile, prewhitening, filter- 
ing, pixelization, and windowing. 

We combine the spectra from individual patches with in- 
verse variance weighting to obtain the final spectrum. 
We also construct a bandpower window function which 
converts a theoretical power spectrum into the corre- 
sponding binned power spectrum which we estimate from 
the maps. We verify these techniques on sky simulations, 
described below. 

4. SIMULATIONS 

Simulations play an important role both in verifying 
that our mapper is working correctly and in understand- 
ing our spectrum. Specifically, we use simulations to test 
for map convergence, estimate bandpower covariances, 
and confirm analytic estimates of the power spectrum 
uncertainties. 

To test for convergence, we perform a set of end-to- 
end simulations of the f 48 GHz and 2f8GHz maps. We 
simulate signal-only time strea ms by mock observa tions 
of input celestial maps (from ISchgal ct al.l I2010bf ) and 
generate noise time-streams which capture the main fea- 
tures of the noise in the data. To test for convergence, 
we simulate a signal-plus-noise map, subtract from it a 
noise-only map at the same iteration and compare the 
power spectrum of the resulting difference map with that 
from a signal-only map. For f 48 GHz, the maps arc 
well converged to multipoles as low as £ > 200 by it- 
eration 1000, and we keep £ > 500 in our final spectrum. 
For the 218 GHz maps, where higher atmospheric noise 
warranted a slightly different mapmaking approach (sec 
Section [2~Tj) . convergence at large scales is slow. After 



fOOO iterations, the 218 GHz maps are converged down 
to £ = 1000, but not at smaller L Moreover, the 218 GHz 
spectrum is highly noise-dominated at low multipoles. 
Thus, we limit spectra involving 218 GHz to i > 1500 in 
our analyses. 

In order to investigate the bandpower covariance and 
to validate the analytic prescription for uncertainties in 
the spectrum, we ran a large set of Monte Carlo simu- 
lations that contain key aspects of the noise properties 
manifested in the data. Due to the iterative nature of our 



Bin index 
20 30 




Fig. 5. — Normalized covariance matrix for decorrelatcd band- 
powers for 148 GHz X 148 GHz spectrum, based on 960 sig- 
nal+noise Monte Carlo simulations. The bins are defined in Ta- 
ble [l] All values on the diagonal are unity by definition. The col- 
orbar has been stretched to reveal the variations in the off-diagonal 
elements. 
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Fig. 6. — Intermediate multipole 500 < £ < 3000 portion of the power spectra plotted with an £ 4 scaling to emphasize the acoustic 
peaks. The thick orange curve s hows the best-fit model including the CMB secondaries and point source contribution taken from our 
companion paper iDunklcv ct al.l 1120101 ). The model depends on the frequency combination. The thin dashed teal line shows the best- fit 
lensed CMB-only theory, which is frequency independent. From top to bottom, the three panels show the 148 GHz, the 148 x218 GHz and 
the 21 8 GHz spectra. For this plot, data spectra and errors from Table [T1 have been scaled by best-fit calibration factors from lDunklcv ct al.l 
J20H]) of 1.022 , 1.02 X 1.09 and 1.09 2 for the 148 GHz, the 148x218 GHz and the 218 GHz spectra, respectively. 



map-making pipeline, it is prohibitively expensive to run 
a large number of end-to-end simulations starting from 
simulated timestreams. Therefore, we adopt an inter- 
mediate approach, in which we find a prescription from 
the data for generating realistic realizations of noise in 
map space, to which we add signal realizations. We start 
with two quarter-season split maps, ?o(n) and Ti(n), 
with corresponding n b s maps labeled n (n) and ni(n). 
We use the difference of these maps to estimate a one-hit 
noise map £>oi(n) (i.e., a noise map representative of the 



variance in each pixel had it been observed only once), 

Dol(ft) = W-r.M (u) 
V n o 1 ( n ) + n i 1 ( n ) 

With the other two splits, we compute another one-hit 
noise map D 2 3(n). We compute the average of the 2D 
power spectra of Dqi and -D23, and use its square root as 
the amplitude for a Gaussian random field, and generate 
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random one-hit noise maps by randomizing the phases 
and converting back to real space. Finally, we divide 
one such map by y/m(ii) to generate a realization of the 
noise for split i. This procedure captures the non- white 
and anisotropic aspects of map noise, as shown by Fig. 01 

We also simulate a random residual point source com- 
ponent in our Monte Carlo realizations. The point source 
model for 148 GHz consists of two components - a radio 
source po pulation and a dusty s ub-mm population, both 
based on iToffolatti ct all (j 19981 ). bourcc maps arc gen- 
erated assuming a Poisson distribution. For each flux 
bin S, spanning the range 0.01 — 20 mJy, we gener- 
ate the number of point sources in that bin as a Pois- 
son realization from the model flux distribution dN/dS 
(the upper limit is chosen to be the approximate 5a 
detection threshold at 148 GHz). Once the number of 
sources in each bin is generated, their individual fluxes 
are dithered to redistribute the fluxes within each bin, 
and the sources are laid down at random positions in 
the map. For 218 GHz we scale the radio and IR com- 
ponents by the ir appropriate spectra l indices, based on 
the findings of iDunklcv ct al.l (|2010l ) . With these set- 
tings, we closely reproduce the level of the Poisson point 
source signal seen in the high multipole regime of the 
data power spectra. We estimate the Poisson contribu- 
tion to the error bars on bandpowers from the scatter of 
500 simulations of point source-only maps, and add that 
as a correction to the (signal+noise)-only error-budget 
obtained from the signal. At the current level of noise 
in our maps, these corrections are largely negligible, and 
we ignore their contribution in the analytic error bar es- 
timation and band power covariance. 

We estimate 148 GHz and 218 GHz spectra for 960 sig- 
nal+noise realizations. Each realization involves simu- 
lating four patches and four splits in each patch for each 
frequency, and estimating the spectra with exactly the 
same methods as used on the data. The resulting spec- 
tra are used for quantifying the covariance between band 
powers. Fig. [5] displays the normalized covariance ma- 
trix, Cav(i,j)/ y / Cov(i, i)Cov(j, j), where i,j denote bin 
indices, showing that correlations between adjacent bins 
are insignificant at the 10% level. 

The uncertainties in the band powers are evaluated us- 
ing the analytic prescription described in the Appendix. 
First, a large set of isotropic white noise simulations are 
generated to validate the analytic formula for the sin- 
gle frequency and the cross-frequency error bars. These 
simulated error bars are found to be in agreement with 
analytic prescription to less than a percent. The next 
set of simulations are run with anisotropic noise, Pois- 
son point source realizations, and with real-space and 
Fourier-space weighting as described in Section [31 and 
small corrections to analytic prescription due to the 
anisotropic nature of the noise and weighting are eval- 
uated against the Monte Carlo simulations. 

5. POWER SPECTRUM 

Applying the methods described in Section [3j we com- 
pute the decorrelated bandpowers for the 148, 218 and 
148 x218 GHz spectra from the unbiased map solutions. 
The bandpowers are presented in Table [1] and are dis- 
played in Figures [S] and [7] In Fig. [5] we have plot- 
ted the lower multipole portion of the power spectrum 
500 < £ < 3000 with an £ A scaling to emphasize the 



higher-order acoustic peaks in the primary CMB spec- 
trum. In the 148 GHz spectrum we clearly resolve the 
second th rough the seve n th pe aks of the CMB. As dis- 
cussed in IDunklcv ct al.l (|2010T ) , the higher-order peaks 
provide leverage on cosmological parameters such as the 
spectral index and its running, the primordial helium 
fraction and the number of relativistic species. Fig. 
on the other hand, emphasizes the high multipole tails 
of the spectra, where the signal is dominated by emis- 
sion from dusty star-forming galaxies and unresolved ra- 
dio sources. The intermediate range 2000 < £ < 4000 
gets a significant contribution from the thermal Sunyaev- 
Zcldovich effect, and lets us constrain the amplitude of 
the SZ spectrum. 

The data are checked for consistency by performing 
various null tests, in which difference maps are con- 
structed to cancel true sky signals, and their power spec- 
tra examined. We describe the suite of tests here. As 
described in Section 13.11 the TODs are split into four 
subsets, generating four maps with independent detec- 
tor noise. We verify that these maps are consistent with 
each other by generating difference maps from each pair, 
and computing the two-way cross spectra from the three 
independent pairs of difference maps: 

c (i-j)x(k-l) = / AT (i-3>(|) AT (fe-0(^ ( 15 ) 

where AT*'"^ = T l — T j and the indices k, I range 
from 0-3. The same point source mask used to mask the 
full datasct is applied to the difference maps before cal- 
culating the spectra; all other settings also remain the 
same. The difference maps, AT^ - ^ , are downweightcd 
by the same n b s maps summed over the four splits used 
in the windowing of the CMB data (see Section 13. 3[) . 
Similarly, when binning the power spectra, we use the 
same azimuthal weights in each patch described in Sec- 
tion 13.31 The three spectra are shown in the top two 
panels of Fig.[8]for the 148 (top left) and 218 (top right) 
GHz maps. We compute the bandpowers in the range 
500 < £ < 10000 (1500 < £ < 10000) for 148(218) GHz. 
The mean of the three spectra is consistent with null with 
X 2 = 42 (14) for 148 (218) GHz, with 51 (31) degrees of 
freedom. 

Another null test probes the consistency of data taken 
with the telescope accelerating as it reverses direction at 
the ends of the scan (turnarounds). Note that for the 
standard maps, no turnaround cut is applied. Four new 
split maps are made cutting data near the turnarounds, 
amounting to losing ~ 24% of the total data. Two differ- 
ence maps are made by pairing one split of the standard 
map with a different split of the new maps (we avoid dif- 
ferencing the same splits as they have very similar noise 
structure), and a two-way cross-power spectrum is pro- 
duced. Any artifact due to the turnaround would be 
left in these difference maps and might produce excess 
power. However, we find the resulting spectra, shown in 
the second row of Fig. [8j to be consistent with null with a 
X 2 = 62 and 18 respectively for the 148, 218 maps (again 
for 51 and 31 degrees of freedom). 

The last two null tests probe possible systematics asso- 
ciated with the focal plane and from diurnal effects. The 
first compares maps made from only the top half of the 
detector array with those from the bottom half. The sec- 
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Fig. 7. — Single and cross-frequency spectra plotted with the conventional 1(1 + 1) scaling, highlighting the behavior at large multipolcs. 
The t hick orange curve sho ws the best-fit model including the CMB secondaries and point source contribution taken from the companion 
paper IDunklcv ct al,l 1 120101 ), The model depends on the frequency combination. The thin dashed teal line shows the best-fit lensed CMB- 
only theory, which is frequency independent. From top to bottom, the three panels show the 148 GHz, the 148x2 18 GHz and the 218 G Hz 
spectra. For this plot, data spectra and errors from Table [T] have been scaled by best-fit calibration factors from IDunklcv ct al.l {201(3) of 
1.02 2 , 1.02 X 1.09 and 1.09 2 for the 148 GHz, the 148x218 GHz and the 218 GHz spectra, respectively. 



ond compares maps processed with data only from the 
middle half of the night to the standard maps. The re- 
sulting spectra are shown in the final two panels of Fig. [8l 
and are also consistent with null. 

6. FOREGROUNDS 

The main foregrounds in the 148 GHz and 218 GHz 
bands at the angular scales considered here are point 
sources. We mask the detected ones as described below. 
Another foreground is the diffuse Galactic dust, which is 
discussed below in so me detail. In the companion paper, 
iDunklev et"aTl (|2010T ). we also consider the contributions 
to the power spectra from the thermal and kinetic SZ 
effects and the clustering of infrared point sources. 

The point sources include radio sources and dusty star- 
forming galaxies (SFG). The radio sources are typically 
active galactic nuclei exhibiting synchrotron-dominated 



spectra, with emission extending down to lower radio fre- 
quencies. The dusty SFG are characterized by absorption 
of ultraviolet photons from star- forming regions by dust. 
The dust reemits into graybody radiation that rises with 
increasing frequency into sub-millimeter bands. 

For the purpose of extracting the underlying CMB sig- 
nal, it is important to identify and mask bright point 
sources in the maps before the power spectrum is com- 
puted. A detailed study of the point source popula- 
tion detec ted in the ACT 148 GHz 2008 survey is pre- 
sented in I MlOl . Point s ources are detected with the 
matched filter method ( Tcgmark fc de O l ivcira-Costal 
119981 : IWright ct all 120091 Tvicira et al.l 120101 ) . To con- 
struct the point source masks, we consider detections 
with SNR > 5 in cither frequency. Over the 296 deg 2 
area used for power spectrum estimation, we find a total 
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Fig. 8. — Null tests of the ACT 148 (left) and 2 18 (r ight) GHz maps. The top row plots the null spectra from difference maps made from 
the four different time splits described in Section 13.11 The second row illustrates the systematic check on whether data near the telescope 
turnarounds arc contaminated, while the bottom two rows check for a systematic gain difference across halves of the focal plane and for 
diurnal effects, respectively. The difference spectra and the chi-squares are computed with the same bins as the standard power spectra. 
For display purposes only, we have re-binned the spectra in this plot with larger bins. The spectra are consistent with null in all cases. 



Predicted Galactic Dust Emission at 219 GHz (fiK) 
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Fig. 9. — Galactic thermal dust emission template at 219 GHz based on IFinkbcincr ct al.l 1119990 "model 8". We cross correlate the 
template patches (numbered here) with ACT patches to look for a correlated dust signal in ACT maps. Patches 0-3 are located off the 
Galactic plane with central Galactic latitudes of —64°, —56°, —43 °, and —28°, respectively. For the dust emission the effective central 
frequency in the ACT 218 GHz band is 220 GHz l|Swctz ct al.ll201Ch . 



of 164 point sources, taking the union of the detections 
at the two frequencies. Before computing the spectrum, 
we mask out a 10' diameter region around each of the 164 
point sources. This amounts to ma skin g about 1% of the 
total area. As discussed in Section [321 the prewhitening 
operation effectively guards against the leakage of power 
due to the application of this mask. 
We estimate the Galactic dust contribution to our 



power spectra by cross-correlating ACT maps with the 
predictions for infrared cirrus emission at CMB frequen- 
cies based on th e mult icomponcnt dust model "8" of 
iFinkbeiner etap (p99h (hereafter FDS). The FDS- 
based dust map s at 148 GHz and 219 GHz are available 
as a part of the iSchgal ct al.l (|2010bl ) simulation^. We 
resamplc these maps into ACT pixclization, and cut out 

These simulations are available at 
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Fig. 10. — Power spectra of dust in the FDS template in the 
patches shown in Fig. pat 150 GHz and 220 GHz in CMB differ- 
ential temperature units. A theoretical lensed CMB power spec- 
trum is shown for comparison (continuous curve) . Also shown is a 
C t <x t~ z curve (dashed) normalized to £ 2 Ci/2tt = 2 (20) liK 2 at 
I = 500 for 150 GHz (220 GHz). Note that a temperature power 
law index of = 1.7 between 150 GHz and 220 GHz corresponds 
to a factor of 10 in the power spectra. 

the ACT patch regions. The 219 GHz FDS map is shown 
in Fig. [SJ Patches through 3 labeled in the figure are 
situated off the Galactic plane with central Galactic lat- 
itudes of -64°, -56°, -43°and -28°, respectively. The 
power spectra of the predicted dust signal in each patch 
are shown in Fig. [TO] 

We proceed with the cross-correlation as follows. We 
express the ACT map as a sum of the cosmic and galactic 
dust components: Tact = 7cmb+C ^fds, where £ is the 
dust amplitude, predicted to be unity by FDS. We take 
two approaches to estimate £: one based on real-space 
operations, and the other based on power spectra. 

In the real space approach, we prewhiten both the ACT 
and FDS maps, convolve the ACT maps with an 6.'1 
FWHM Gaussian beam appropriate to the FDS resolu- 
tion, apply the ACT point source mask described above, 
filter out modes below I = 500 (1500) for 150 GHz (220 
GHz) and produce the maps Tl CT (n) and 7j? DS (n) in 
patch i. Then we estimate £ in patch i as 

Jd 2 nri CT (n)7g DS (n) 

The error on the above estimate is obtained by cross- 
correlating 960 random CMB plus noise simulations (see 
Section with the dust maps. In the spectrum-based 
approach, we estimate £ in each multipolc bin b as the ra- 
tio of the ACT-FDS cross-spectrum to the auto-spectrum 

http : //lambda. gsf c .nasa. gov/toolbox/tb_cmbsim_ov . cf m 

To be precise, we use the actual frequencies tor which the FDS 
maps were made, and the effective dust frequency for ACT, but in 
practice, the differences are negligible. 



of the FDS map for the same patch: 

£,ACTxFDS;i 
£b = ,-,FDSxFDS;i " 

The ACT-FDS cross-spectrum is calculated using 
both pixel-space weighting and Fourier-space azimuthal 
weighting, as with the ACT power spectrum, with the 
FDS and ACT beams deconvolved. The final estimate, 
£s pcc is obtained as the mean across the bins over the 
range 500 < b < 3000 (1500 < b < 3000) for 150 GHz 
(220 GHz) where the upper limit is dictated by the FDS 
resolution. The uncertainty on £g pec is estimated as the 
scatter in the same quantity computed from the cross- 
spectra between the FDS maps and 960 random CMB 
plus noise realizations. As a check for systematics, we 
circularly shift the FDS patches by one patch, such that 
patchJ) of ACT is crossed with patch_3 of FDS, patch_l 
of ACT with patchJD of FDS etc., and recompute the real 
and Fourier space estimates of £. The results are shown 
in Fig. 1111 The errors shown arise from treating the FDS 
template as deterministic. 

This figure shows that a dust-correlated signal is ob- 
served. More work will be needed to understand how well 
the FDS template r epresents faint dust we ll off the Galac- 
tic plane (see, e.g.. IVencziani ct al.ll2010T ). The fact that 
such a small signal (£ 2 Ci/2n ~ 0.1 - 1.0 /iK 2 at I = 3000 
and 148 GHz) can be recovered with our map-based 
technique demonstrates the power of our max i mum - 
likelihood map estimation. Recently lHall et ail (|2010D 
have also reported a measurement of infrared cirrus in 
the SPT maps. 

Subtracting the FDS template from our maps has neg- 
ligible effect on the power spectra, consistent with the 
minuteness of the dust signal. 

7. LENSING OF THE CMB 

The CMB photons are deflected by large-scale struc- 
ture potentials along their pat h from the last scattering 
surface at z ~ 1100 to us (see iLewis fc Challinod[2006L 
for a review) . The typical (rms) deflection in the ACDM 
model is about 2 '7 and the deflection pattern is coherent 
over degree scales, comparable to the size of the acous- 
tic features on the primary CMB. These coherent deflec- 
tions produce distortions of the hot and cold spots on 
the CMB, leading to a broadening of their size distribu- 
tion. In the power spectrum, this effect manifests itself 
as the smoothing of the acoustic peaks, which can be 
used as a signal to look for lensing. The first attempt at 
a detection of lensing in the power spectrum was made in 
iReichardt ct all (|2009bl ) using the data from the ACBAR 
experiment. They quantified the effect of including lens- 
ing in their analysis through the log-ratio of the lensed 
to unlensed Bayesian evidence (Aln£), and by combin- 
ing WMAP5 an d ACB AR datasets found Alnf = 2.63. 
iCalabrese et all (|2008D analyzed the ACBAR data with 
a different approach, where they introduced a scaling 
of the power spectrum of the lensing potential (Cf — > 
A^Cf , with Al = corresponding to no lensing and 
Al = 1 to the standard ACDM expectation). With 
WMAP5+ACBAR they found A L = 3.0 ±g;| (68% confi- 
dence le vel) . According to footnote 17 of lRcichardt ct all 
(l2009b[ ). a reanalysis of the ACBAR data using this pa- 
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Fig. 11. — Left: Amplitude, £, of dust emission in ACT patches relative to FDS predictions for 150 GHz (open symbols) and 220 GHz 
(closed) symbols. Results obtained through the real-space and Fourier-space techniques are shown with circles and diamonds, respectively. 
Right: Results of a null test performed by shuffling around the FDS patches, so that patch of ACT is crossed with patch 3 of FDS etc. 
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Fig. 12. — Lensing of the CMB smooths out the acoustic peaks 
in the CMB power spectrum. The best fit model with lensed CMB, 
secondaries, and point sources is shown as the thick orange curve, 
while the same with no lensing is shown with the thin green curve. 



Our 148 GHz power spectrum is shown against lensed 
and un lcnscd models in Fig.lTSl We use the parametriza- 
tion of ICalabrcsc "eFaLl (pOOa T and using the parameter 
esti mation method o logy d escribed in the companion pa- 
per iDunklcv et al.l (|2010D , we constrain the lensing pa- 
rameter Al based on WMAP7 and ACT power spec- 
tra. Fig. [13] shows the marginalized ID likelihood for 



.4 



1.3 



+0.5(+1.2) 



L using WMAP7+ACT. We hnd Al - 

at 68% (95%) confidence, with the best-fit lensed CMB 
spectrum with Al = 1 having an improved goodness-of- 
fit to the WMAP7+ACT data of A x 2 = 8 less than the 
unlenscd model. 

We check for systematics that might have given rise 
to a spurious lensing signal. The projection scheme 
(cylindrical-equal-area) used for the ACT maps is not 
particularly optimized for lensing studies — so we test 
whether this projection could introduce a lensing-likc 
signal. We simulate a low-noise unlensed CMB sig- 



nal and run it through the mapmaking pipeline, and 
try to reconstruct a "lensing convergence" in the re- 
sulting map using standard quadratic estimator tech- 
niques (|Hu fc Okamoto|[2002t ) . We find the reconstructed 
convergence power spectrum to be consistent with null 
(there is a small known bias at high multipoles that we 
have entirely traced to mode-coupling due to the finite- 
ness of the patch), showing that the projection does 
not introduce any significant lensing-like signal in our 
patches. To further test if any other step in our pipeline 
could produce spurious peak smearing in the power spec- 
trum, we generate an end-to-end simulation of a noisy 
map exactly as described in Section 21 only this time re- 
placing the lensed CMB signal time-stream with its un- 
lcnscd version. The resulting maps are then processed 
through the power spectrum pipeline, and the power 
spectrum is then analyz ed with the pa r amete r estima- 
tion method described in IDunklcv ct all (|2010f ) . We find 
the lensing amplitude parameter, Al, described above to 
be consistent with zero (see Fig. [T3J) • 

Here we have used the smearing of the acoustic features 
to look for the lensing signal. More promising ways of 
extracting this signal involve using information beyond 
the two-point correlation (optimal quadratic estimators) 
and cross-correlation with tracers of large-scale structure. 
Early efforts cross-correlating the lensing reconstruction 
in WMAP 3-year maps with luminous red galaxies, ra- 
dio s ources and quasars (|Smith ct al.H2007t Iffirata ct all 
12008(1 yielded ~ 3c evidence for lensing. Data from 
WMAP are not the best suited for this kind of study be- 
cause of its limited angular resolution. With arcminute 
resolution CMB data, such as from ACT, there is a much 
higher potential for a detection. Such efforts are under- 
way with the ACT maps. 

8. DISCUSSION 

We have presented a measurement of the power spec- 
trum of the CMB sky observed with the Atacama Cos- 
mology Telescope at 148 GHz and 218 GHz. The 148 GHz 
spectrum spans a large dynamic range from I = 500 to 
£ = 10000, covering the damping tail, where the primary 
anisotropics with the higher-order acoustic peaks domi- 
nate, to the composite high multipolc tail of the CMB 
where emission from dusty galaxies and radio sources and 
the SZ effect contribute. The second through the seventh 
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TABLE 1 

Single frequency bandpowers 
B b =l b (t b + l)C b /2n (nK 2 ) 



148 GHz 148 GHz X 218 GHz 218 GHz 



t range 


cent ral 


Ob 




o b 


rr(R, \ 

&\pb) 


o b 


rr(R, \ 


491 - 


540 


515 


2423.0 


270.8 










541 - 


590 


565 


2701.4 


274.8 










591 - 


640 


615 


1952.1 


205.4 










641 - 


690 


665 


1701.9 


167.7 










691 - 


740 


715 


1895.7 


179.5 










741 - 


790 


765 


2386.5 


202.0 










791 - 


840 


815 


2415.5 


200.5 










841 - 


890 


865 


2168.0 


178.2 










891 - 


940 


915 


1567.5 


138.3 










941 - 


990 


965 


1353.6 


109.9 










991 - 


1040 


1015 


976.0 


89.2 










1041 - 


1090 


1065 


1089.5 


92.9 




_ 






1091 - 


1140 


1115 


1229.6 


100.2 










1141 - 


1190 


1165 


1169.7 


88.2 










1191 - 


1240 


1215 


861.9 


73.4 










1241 - 


1290 


1265 


819.3 


68.7 










1291 - 


1340 


1315 


740.8 


60.6 










1341 - 


1390 


1365 


746.8 


60.5 










1391 - 


1440 


1415 


782.4 


64.4 










1441 - 


1490 


1465 


718.8 


55.8 










1491 - 


1540 


1515 


619.3 


50.1 


559.6 


43.6 


585.9 


87.6 


1541 - 


1590 


1565 


495.6 


44.5 


489.0 


40.2 


508.0 


83.0 


1591 - 


1640 


1615 


467.4 


38.8 


428.3 


35.3 


480.1 


77.9 


1641 - 


1690 


1665 


380.4 


34.0 


382.0 


32.7 


446.5 


74.0 


1691 - 


1740 


1715 


351.6 


33.6 


347.8 


32.1 


436.8 


75.0 


1741 - 


1790 


1765 


359.6 


32.3 


336.7 


30.7 


470.6 


74.1 


1791 - 


1840 


1815 


334.4 


29.4 


336.9 


28.3 


396.1 


65.3 


1841 - 


1890 


1865 


274.8 


25.3 


276.6 


25.8 


371.5 


65.0 


1891 - 


1940 


1915 


226.3 


23.3 


198.8 


23.7 


199.7 


62.1 


1941 - 


1990 


1965 


259.7 


23.1 


254.6 


23.9 


310.2 


61.9 


1991 - 


2040 


2015 


235.8 


21.6 


226.0 


22.5 


340.2 


60.6 


2041 - 


2090 


2065 


242.8 


22.8 


227.7 


23.9 


283.9 


63.5 


2091 - 


2140 


2115 


168.8 


18.6 


206.9 


21.2 


268.3 


57.8 


2141 - 


2190 


2165 


164.7 


17.8 


168.4 


19.7 


203.4 


54.9 


2191 - 


2240 


2215 


152.2 


17.2 


151.0 


19.7 


269.7 


60.5 


2241 - 


2290 


2265 


132.2 


15.8 


114.6 


18.1 


206.7 


56.3 


2291 - 


2340 


2315 


114.9 


14.7 


136.5 


18.1 


135.8 


54.2 


2341 - 


2390 


2365 


106.3 


15.4 


135.3 


19.1 


222.3 


59.2 


2391 - 


2700 


2550 


82.9 


5.4 


99.7 


6.7 


157.9 


21.2 


2701 - 


3000 


2850 


57.5 


5.3 


72.6 


6.5 


189.1 


20.8 


3001 - 


3400 


3200 


37.2 


4.8 


54.0 


5.5 


158.3 


16.9 


3401 - 


3800 


3600 


39.8 


5.6 


50.9 


5.9 


145.1 


16.6 


3801 - 


4200 


4000 


48.1 


6.8 


62.9 


6.7 


204.5 


17.9 


4201 - 


4600 


4400 


34.9 


8.3 


74.6 


7.8 


198.7 


19.7 


4601 - 


5000 


4800 


36.4 


10.2 


98.0 


9.3 


248.7 


22.2 


5001 - 


5900 


5450 


55.8 


9.6 


106.3 


8.3 


281.7 


18.9 


5901 - 


6800 


6350 


56.6 


15.4 


131.7 


12.5 


378.3 


26.8 


6801 - 


7700 


7250 


127.8 


25.6 


127.7 


19.1 


463.7 


38.1 


7701 - 


8600 


8150 


115.3 


42.4 


228.0 


30.0 


620.6 


55.6 


8601 - 


9500 


9050 


208.3 


73.1 


175.9 


48.6 


812.8 


85.0 


9501 - 


9900 


9750 


316.9 


172.4 


152.4 


108.1 


695.0 


176.1 



acoustic peaks of the CMB are clearly visible in this spec- 
trum. For the 148 x218 GHz and 218 x 218 GHz spectra, 
we present measurements from 1500 < £ < 10000. The 
high multipole information from these spectra help us 
constrain the nature of the point source populations that 
contribute to this range, and in turn, constrain the SZ 
contribution to the 148 GHz spectrum. We measure the 
presence of a very faint dust signal at high Galactic lati- 
tudes. Recovering such a faint signal gives us confidence 
in the fidelity of the maps. We find evidence of gravita- 
tional lensing of the CMB in the power spectrum at the 
2.8tr level. Constraints on cosmological parameters and 
the interpretation of high multipole spectra in terms of 
point source populations a nd the SZ effect are pr esented 
in the companion paper bv lDunklev claTl (pOlOl ). 
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APPENDIX 



CROSS-FREQUENCY CROSS-SPECTRUM ERRORBARS 



Here we derive an analytic expression for the expected error bars on the cross-frequency cross-power spectrum, 
assuming the signal and the noise can be approximated as Gaussian random fields. We do not treat the Poisson 
component from point sources here, as the Poisson noise correction is negligible at our current noise level. As before, 
we denote the two frequency channels as A and B and use lowercase romans k, I, etc. to denote the sub-season data 
splits. 

The estimator for the cross-frequency cross-power spectrum in bin b is given by: 



AxB 



1 



nd{nd - 1) 



■y^ Q(iAxjB) 



where rid is the number of sub-season data splits (four in our case), and 



C, 



(iAxjB) 



= -E T * 



ArpjB 



(Al) 



(A2) 



where the sum is over pixels in the annulus in Fourier-space contained in the bin b, and v\, is the effective number of 
pixels in the annulus over which the sum is taken. Evidently, 



(iAxjB) 



AxB 



1 Cb, 

c b . 



(A3) 
(A4) 



First, we will compute the expected covariance between a pair of estimators: 



a (iAxjB);(kAxlB) 
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(iAxjB) 



(kAxlB) 



C, 



(kAxlB) 
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where i ^ j and k ^ I. From the definition (|A2[) . we can write 



e 



(iAxjB);(kAxlB) 
bb 



1 



EE[ 



rp*lArpjB rp*kArpll 
J- f J-g J-£r -L pi 



b lEbl'eb 

Applying Wick's Theorem to the above formula, we get, 



£,(iAxj\B)\ /^(kAxlB) 
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bb 



ry (iAx]B)\ / ^j(kAxlB) 



v b\C b 



+ u b (c { ^ kA) ) (c^ lB) ) + v b (c^ AB) ) (C, 



(kAxjB) 



(iAxjB) 



C, 



(kAxlB) 



1 r 

Vb 



c (iAxkA)\^ / c QjBxlB)\ + ^ c (iAxlB)^ ^(kAxjl 



(A6) 
(A7) 



Now we will examine various cases of the above variance. 
Type 1: All four splits different, i ^ j ^ k ^ I. In this case, all four ensemble averages are cross-power spectra and 
evaluate to the underlying spectrum C b , giving 



n (iAxjB);(kAxlB);Type 1 
^bb 



Vl, 



r 2 



(A8) 



Type 2a: One split in common, of channel A, (i = kj^l^j). Here, the first ensemble average in (|A7[) is an 
auto- frequency auto-power spectrum which evaluates out to (Cb + N AA ), where N AA is the noise in the channel A 
data split. This gives 



e 



(iAxjB);(kAxlB);Type 2a 
bb 



Vl, 



-{2C£ + C b N{, 



AA\ 



Type 2b: One split in common, of channel B, (I = j ^ k ^ I 



^(iAxjB ; (kAxlB ;Type 2b 1 {or ,2 , n Ar BB\ 
V bb = — (2C fc + C b l\ b ). 

Vb 



Type 2c: One split in common, one of channel A, the other of channel B (i = I or k = j ). 

Q gAx 3 B);(kAxlB);Typc 2c = + CbN AB^ 



Vl, 



where N AB represents any noise that is correlated across the two frequency channels. 
Type 3a: Two splits in common [i = k, j = I ). In this case, we have 



e 



(iAxjB);(kAxlB);Typc 3a 
bb 



1 

Vb 



2C z b + C b (N AA + N^) + N b AA N£ 



BB\ 



T AA atBB-I 



Type 3b: Two splits in common (i = I, j = k). This evaluates to 

e (,Ax,B);(/cAx iB );Typc 3b = J_ ^2 + 2Cb N AB + N AB N AB ] 



Vl, 



Next, we turn our attention to the variance of the mean spectrum estimator, C^ xB , defined in ([All) : 



0mean;mean ^ | ^j/i x n / 
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Expanding one of the terms out, we have 

0mcan;mcan \ ^ / / 
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(All) 



(A12) 
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(A14) 



(A15) 



Each of the n d (n d — 1) terms in the sum is statistically equivalent to the other and should evaluate to be the same. 
Hence, 

0mcan;mcan Q(iAx j.B);mean (A16) 



which can be expanded out as 



e 



mcan;mcan 



E e 

d-D .4?.. 
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nd(n d - 1) 
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The last piece of the calculation is to figure out how many times each of the types defined above appear in this sum. 
Since i and j are fixed, it is apparent that for Type 1, k can be chosen in (rid — 2) ways, and I in (rid ~ 3) ways. So, 
Type 1 appears (rid — 2)(rid — 3) times. For Type 2a k is fixed by equality with z, so I can be chosen only in (rid — 2) 
ways. Therefore, Type 2a appears (rid — 2) times. From symmetry, Type 2b also appears (rid — 2) times. For Type 2c, 
either i = I or k = j which gives 2(rid — 2) terms. Finally, for each of Type 3a and 3b, all four indices are completely 
fixed, so they appear only once. Putting these together, we have: 
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which simplifies to 

Qincan;mean 

Note that if there is only one channel A 
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B, this reduces to the familiar form (see Eq. (9) in lFowlcr ct al.ll2010f ) 
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Finally, if we have n p patches with equal noise, the above variance should be divided by n p . In practice, a weighted 
combination of the individual patch variances is used. Under the assumption of isotropic noise and filtering, the effect 
of the data window, W , can be taken into account by correcting the number of modes as v b v b w\jw^ where w n 
represents the n-th moment of the window. There is a small additional correction due the anisotropic nature of the 
noise and filtering, which we calibrate against Monte Carlo simulations 

For estimating parameters from the three spectra, C\ 
written as 



148x148 ^148x218 j ^218x218 ru i r ,. 

C b and C b a joint likelihood function is 
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where 
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,A-s,c-n 



= a 



^AxB 



with 



and can be computed in a similar fashion as above. 
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